Analysis of the loss of coherence in interferometry with macromolecules 
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We provide a self-contained quantum description of the interference produced by macromolecules 
diffracted by a grating, with particular reference to fullerene interferometry experiments 0-0 We 
analyze the processes inducing loss of coherence consisting in beam preparation (collimation setup 
and thermal spread of the wavelengths of the macromolecules) and in environmental disturbances. 
The results show a good agreement with experimental data and highlight the analogy with optics. 
Our analysis gives some hints for planning future experiments. 
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I. INTRODUCTION 

The aim of this paper is to provide a theoretical deriva- 
tion of the measured beam intensity profile for interfer- 
ometry with macromolecular beams, under the influence 
of the processes inducing loss of coherence, consisting 
both in beam preparation (collimation setup and ther- 
mal spread of the wavelengths of the macromolecules) 
and in environmental disturbances. 

We have been motivated by the impressive experiments 
with fullerene made by Zeilinger's group 0,0- I n these 
experiments, thermally produced beams of heavy macro- 
molecules are collimated, diffracted by a grating and then 
detected on a distant screen. The diffraction pattern so 
produced shows the typical interference profile of wave 
phenomena in the presence of incoherent contributions, 
and the reduced fringe visibility observed reminds very 
much of Kirchhoff diffraction with thermal light produced 
by an extended source. By taking into account the ef- 
fects of the interaction of the macromolecules with the 
environment and with the photons they emit by inter- 
nal cooling, we shall provide a self-contained quantum 
description of these experiments, that does not rely on 
methods of classical optics. 

Our analysis is based on two main ingredients. The 
first one is a matter of principle: it is the formula for the 
statistics of particle arrival position and time on a distant 
surface given by Eq. below. According to this for- 
mula, the intensity pattern revealed on a distant screen 
can be expressed in terms of the large-distance asymp- 
totic behavior of the time-integrated quantum current 
(see Eq. © below). This formula was conjectured long 
time ago within the framework of scattering theory Q, 
and only recently has been derived 0, 0, El i extended to 
the mesoscopic regime, and physically motivated within 
the framework of Bohmian mechanics |7( . 
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The second ingredient is a matter of analysis to sim- 
plify the dynamics of a test particle moving in a quan- 
tum medium: it is the model of Joos and Zeh y| for the 
phenomenological description of processes inducing loss 
of coherence in quantum systems. In this model, the re- 
duced density matrix of the system evolves autonomously 
according to a "Boltzmann-type" master equation. The 
effect of the environment is summarized by "a collision 
term" , added to the free dynamics of the system, which 
takes into account the decoherence, i.e., damping of the 
off-diagonal terms of the density matrix in position rep- 
resentation. 

Starting from these premises, and by means of some 
approximations that are reasonable in the common ex- 
perimental conditions for interferometry with heavy par- 
ticles, we shall derive an easy relation useful to describe 
diffraction patterns. 

Finally, we shall provide a theoretical fit for the exper- 
imental data reported in 2], and we shall discuss some 
predictions about the dependence of the interference pat- 
terns on physical parameters such as the mass of the 
molecules, the pressure at which the experiment is per- 
formed and the distance of the detection screen. These 
predictions can be of some relevance for planning new 
experiments. 



II. MEASURED INTENSITY 

In the experiments with fullerenes a thermally pro- 
duced beam of heavy macromolecules is collimated, 
diffracted by a grating and then detected on a distant 
screen (see Fig. QJ. The grating is composed by parallel 
slits, and it is periodic of period D. The grating and the 
detection screen lie in parallel planes, which are orthogo- 
nal to the longitudinal direction of beam propagation (y 
direction), the so-called optical axis. The screen is placed 
at the distance L from the grating. During the flight from 
the grating to the screen the fullerenes interact with air 
at low pressure, as well as with thermal photons, and get 
entangled with the photons emitted by relaxation of the 
internal excited states. What is the theoretical prediction 
for the particle intensity measured at the screen? In order 
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to answer to this question we shall proceed in two steps: 
firstly, in Section fll Al we shall consider the case in which 
the dynamics of center of mass of the diffracted particle 
is governed by free Schrodinger evolution. Though this 
approximation is unrealistic for fullcrencs, it is a good 
approximation for lighter particles as electrons or neu- 
trons. Then, starting with Section Til Bl we shall refine 
the description and take into account the influence of the 
environment on the motion of the particle. 



A. Free evolution 

In a typical detection experiment the count statistics is 
obtained by summing a large number of events in which 
the particle crosses the detection screen at a random 
time. What is the appropriate quantum prediction for 
such a statistics? This question is not quite as innocent 
as it sounds; it concerns in fact one of the most debated 
problems in quantum theory, the problem of time mea- 
surement, specifically the problem of arrival time, and 
position at such time. It is well known that there is no 
self-adjoint time observable of any sort, and since the ar- 
rival position is the position of the particle at a random 
time, it cannot be expressed as a Heisenberg position 
operator in any obvious way. Bohmian mechanics does 
provide, however, a remarkably simple answer (for an up- 
dated review of Bohmian mechanics see Q and references 
therein): Let S be a surface in physical space, \&(r, i) be 
the wave function of a particle, and 

J(r,t) = Ai m [*(r,t)* V r *(r,<)] 

be the associated quantum current satisfying the conti- 
nuity equation 



Detection screen 



d\*(r,t)\< 
dt 



+ V r • J(r,t) = 0. 



Then the joint probability Prob(R T 6 dS, T <E dt) that 
the particle crosses the surface element dS of the surface 
S at the point r in the time between t and t + dt is given 
by 



Prob(R T £ dS, T 6 dt) = J(r, t) ■ dS dt 



(1) 



provided the current positivity condition J(r, t) ■ dS > 
is satisfied (a condition on both the wave function ^> and 
on the surface S) . See Q for a general derivation of this 
(for applications to mesoscopic physics see 0, 0, ^| , 
in this regard, see also fl^l. 

We would like now to apply such a probabilistic pre- 
diction to the typical situation of a diffraction experi- 
ment with particles of mass M diffracted by a grating 
and then detected on a distant screen. The geometry is 
that of Fig.^and, by (approximate) translational invari- 
ance along the slit extension — the z axis — without loss of 
generality, we can consider the dynamics to be effectively 
two dimensional: r = (x,y), where x is the coordinate 
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FIG. 1: Geometric configuration of the diffraction grating and 
the detection screen. 



along the grating perpendicular to the slit axes and y is 
the coordinate perpendicular to both the grating and the 
detection screen, i.e., along the so-called optical axis. 

Let us make the physically reasonable assumption that 
the initial (t = 0) wave function produced by the grating 
factorizes 

*(r,0) = ^(:r) My), 

with the size Ax of the support of tpo being that of the 
grating. Then \1/ evolves freely according to Schrodinger 
equation 



9* 



2M 



Vr 2 * 



(2) 



until the particle is detected on the screen in the xz- 
plane placed at a "large" distance y = L. This last 
condition — see below for a suitable specification of how 
large L should be — ensures that the current positivity 
condition is fulfilled. Then, according to Eq. QJ, the 
probability density that the particle crosses the screen at 
the point x is 



f + OO 

I(x) = I dt J y (x,y,t) 



(3) 



y=L 



where J y is the longitudinal component of the quantum 



I(x) 



h 

M 



dt \ip{x,i)\ l Im 



dy 



y=L 



(4) 

For a large ensemble (beam) of particles identically 
prepared in the same initial state, by the law of large 
numbers, I(x) is proportional to the local intensity mea- 
sured at the screen — and without loss of generality, the 
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proportionality constant, which is easily determinable by 
the total count statistics, will be hereafter assumed to 
be 1. 

Let us now make some simplifications: Assume that 
the momentum component p y is sharply defined, i.e., 
Ap y <C p y , so that with the initial wave function is asso- 
ciated a well-defined de Broglie wavelength 



A - h/p y < Ay . 



(5) 



Then we may approximate Schrodinger's evolution of 4>o 
with a classical propagation at velocity v — p y /M, so 
that I(x) gets approximated by 



I\{x) — v I dt \ip(x, t)\ I (f>o (L — vt) | 



(6) 



Suppose furthermore that the detector distance L is 
much larger than the position spread in the longitudinal 
direction, 



L > Ay. 



(7) 



Then the time-integration in Eq. 10 gives appreciable 
contributions only for 



L MX r 

t = T = - = Ti; L - 

v Zirn 



(8) 



which is the so-called "time of flight," that is the time 
spent by the particle to reach the detector. Thus, Eq. © 
can be further approximated as 



h(x) = U(x,T)\ 



(9) 



where ip(x,t) is the solution of one-dimensional free 
Schrodinger's equation, i.e., 



tp(x,t) = 



M 



2niht 

Consider now \ip(x,t)\ 2 , i.e., 
M 



dx e^ {x ~ Xo> ipo(xo) . 



2-Kht 



dx dx' e miz 2 o-x'o+^(*'«-x«)] ip ( Xo ) ip ( x ' y 



and note that in the integrations both xo and x' are 
bounded by Ax, the support of i/jq. Thus, if 



Af(Ax) 2 
M 



< 1 



(10) 



we have that emt fao x 'o. « 1 and therefore that \ip(x, i)| 2 
is approximated by 



dx dx' e^o-^VoOo) M4Y 



M 



One can easily recognize that the above expression (mod- 
ulo an overall proportionality constant) is nothing but 
the square of ^{k), the Fourier transform of V'o(x), com- 
puted for k — Mx/(M). Therefore, by replacing in 



Eq. © such an approximation of |?/>(2;,i)| 2 for t = T, 
and recalling definition JSJ for T, we arrive at a further 
approximation for the intensity I(x), namely, 



h{x) = XL 



ipo 



f2nx\ 
{-XLj 



(11) 



It is important to observe that the regime in which this 
approximation holds is that fixed by Eq. i|10|) for t = T 
i.e., the regime characterized by the condition 



Ax X 
T ^ Ax" 



(12) 



Note that condition l|12|l is indeed Fraunhofer's condition 
of classical optics and Eq. (|llf> is the corresponding for- 
mula for the intensity of classical Fraunhofer diffraction 
theory [l4|. according to which the large- distance inten- 
sity of a diffracted field is the squared modulus of the 
Fourier transform of the field distribution on the diffrac- 
tive grating. (For example, for a double-slit diffractive 
grating, with aperture d of each slit, distance D between 
the slits, and with ipo being the characteristic function of 
the slits, Eq. (fill) becomes the standard formula for the 
intensity in the double-slits experiment, namely, 



1 o . 2 
-smc 



ndx 
~XL 



f2irDa 



V XL 



where Iq is the intensity detected for x = and, as usual, 
sinc(x) = x^ 1 sin(x).) 

In this regard, it should be observed that condition 
(|Trj|l leading to ^(x) is a known condition nm. it 
corresponds to a "large time" regime 



T> T 



M(Ax) 2 



Ax 



Ax 



Ap x /M Av x 



(13) 



for which the solution of Schrodinger's equation is ap- 
proximated by 



ip(x,t) 




-(f) 



Such a wave is what in 15] has been named "local plane 
wave": a wave that locally looks like a packet, having 
amplitude and local wave number that are slowly varying 
over distances of the order of the local de Broglie wave- 
length. Such waves are associated with classical motion 
of particles and a rough estimate of the time needed for 
the formation of such waves is indeed the time r above 
[l5| . So, evidence to the contrary notwithstanding, the 
particle motion in the Fraunhofer region, that is, the par- 
ticle motion on the time scale l|13[l . is indeed classical 
motion. 

It might be useful to compare the domain of validity of 
the various approximations. Approximation i2(x) holds 
under the spatial condition Q , that is on the time scale 

MX , 
T»— Ay. 
2nn 
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The temporal and spatial conditions leading to the 
Fraunhofer-like approximation I§(x) are deduced by 
Eq. l|lQp for t = T. For an initial packet with 
Aa; ~ Ay, comparison of Ax <C L (leading to l2(x)) with 
Aa: <C V XL (leading to h(x)) shows that the Fraunhofcr 
approximation is realized on much larger space and time 
scales. 

A final remark: It could be objected that it is indeed 
© the basic formula for the statistical predictions of 
detection experiments — after all, it is this formula that 
seems to correspond directly to the standard statistical 
interpretation of the wave function. This objection, how- 
ever, misses the point altogether: I-2,{x) is only an ap- 
proximation; the time at which the particle crosses the 
screen is typically random, and it can be treated as the 
deterministic quantity given by the time of flight only 
when condition J7J) is satisfied. In a sense, it is true that 
in the regime of large distances there is no experimen- 
tal difference between I(x) and the approximations we 
have considered (indeed, this is a consequence of what it 
has been proven with great generality in 0, 0)- How- 
ever, experimental research on near field interferometry 
may explore regimes in which these approximations fail, 
e.g., when statistical fluctuations in the arrival time be- 
come experimentally relevant, as we shall comment in 
Section IVI CI Hence the need of an exact formula for 
the intensity. And, while the standard quantum formal- 
ism fails to provide the exact expression for the intensity, 
the formula given by Eq. , which clearly looks right, is 
naturally deduced from first principles of Bohmian me- 
chanics 0. 



B. Interaction with the environment 

In the more general case of a quantum particle inter- 
acting with its environment, the evolution cannot be any 
more treated in terms of one-particle Schrodinger equa- 
tion because entanglement can be fastly developed. In 
this case, the statistical predictions concerning experi- 
ments performed on the particle are governed by the re- 
duced density matrix p(r, r'), which is obtained from the 
wave function describing the particle and its environment 
by integrating out the configurational degrees of freedom 
of the environment; for example, for an environment with 
N particles and total wave function ^(r, ri, . . . , r/v), the 
reduced density matrix is given by 

p(r, r') = /dri . . . drjv^r, r X) . . . , r N ) *(r', r x , . . . , r N )*. 

(14) 

As far as detection experiments are concerned, the fol- 
lowing natural generalization of Eq. can be put for- 
ward: The joint probability that the particle crosses the 
surface element dS of the surface S at the point r in the 
time between t and t + dt is still given by Eq. , but 



with current J given now by 



J(r,f) = ^Im [V r p(r,r',i)] r , = 



(15) 



A detailed derivation of this result will be given elsewhere 
|H^ . In this regard, a key observation is that J given by 
Eq. I|15fl is indeed the right probability current entering 
in the continuity equation for the probability density of 
position p(r, r, t), 



dp{r,r,t) 
dt 



V r • J(r,t) = 0. 



We may now make more realistic the analysis of Sec- 
tion lll Al bv allowing that, during the flight from the grat- 
ing to the screen, the particle of mass M diffracted by 
the grating interacts with particles of the environment 
(say, air molecules). As before, the geometry is that of 
Fig. ^ so that the dynamics is effectively two dimen- 
sional, i.e., as before, r = (x, y). Though the initial state 
produced by the grating needs not be anymore a pure 
state, we maintain the physically reasonable assumption 
of factorization at t = 0, 



p{r,r',Q) = p^\x,x')p^{y,y'). 

For an environment of N particles, the time evolu- 
tion of p(r,r',t) is that induced, according to Eq. lfH|l . 
by the Schrodinger evolution of total wave function 
# = *(r,r u ...,r N ,t), 

f)\p h 2 

iH ~dt = ~2M Vr + + Hiat *' (16) 

where Hg nv is the total Hamiltonian of the N particles 
(the sum of the kinetic and potential energies), and H mt 
is the interaction potential between the particle and the 
other N particles. 

Accordingly, the exact formula for the probability den- 
sity that the particle crosses the screen at the point x 
is given by Eq. (|3J), where J y is now the longitudinal 
component of the probability current given by Eq. I|15|) . 
We shall now simplify the expression for the intensity, in 
analogy with the treatment of Section III Al by exploit- 
ing the typical physical conditions of interference exper- 
iments. 

Motion along the y direction is typically "very fast," 
being characterized by a "very short" wavelength A, much 
smaller than all the other relevant lengths scales (such as 
the spreads Ax, Ay, and the screen distance L). Accord- 
ingly, we have an effective preservation of the factoriza- 
tion of the initial state, and Eq. becomes 



(y), 



h f +oa 



y=L 



Note that, due to the condition of fast motion along 
y direction, we may assume wave-packet motion, i.e., 
p^ v \y, y', t) = (j)(y, t) 4>(y' , t)* , and consider the evolution 
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of the wave packet <j> to be classical. Thus, proceeding 
as in Section Hi Al in going from Eq. @ to Eq. JJjJ, and 
a part from a caveat we shall discuss below, we arrive at 
the following approximation for the measured intensity 

I 2 (x)=pW(x,x,T), (17) 

where T is, as before, the time of flight given by Eq. JBJ). 

The caveat is the following: in the free case a crucial 
condition for the validity of Eq. © is that A <C Ay. In 
case of environmental interaction the momentum spread 
increases due to scattering events and the previous con- 
dition is no more sufficient to assure a well defined de 
Broglie wavelength along the longitudinal direction. An 
important consequence of the analysis of Section is 
that interaction with the environment produces an ef- 
fective reduction of the relevant length scales over which 
quantum coherence is preserved and this reduction is con- 
trolled by what we shall call the "coherence length" and 
denote by I. In general, the validity of Eq. l(T7|) is assured 
by A <§C Ay and A <C £; for relevant incoherence effects, 
we have t < Ay and thus the crucial condition becomes: 
A < I. 

Let us consider the case of fullerene, and let the initial 
state be the state at the moment of the splitting produced 
by the diffractive grating. It turns out that, at this time, 
the motion of fullerene along the y-direction can be de- 
scribed by a narrow wave packet translating with veloc- 
ity v. In fact, according to Tab.[IJ the typical de Broglie 
wavelength for fullerene is A w 10~ 12 m. The analysis 
performed in Section IIII Al leads to I m 10~ 7 m (see 
Tab. IHIfl . whence A <C £ (and since L « lm, Ay <C L). 
Thus, Eq. (|17|l provides a good approximation for the 
measured intensity of fullerenes. 

N.B. Eq. (jl7|l is the basic equation of this paper. In 
order to avoid notational complexity, when no confusion 
will arise and unless otherwise stated, we shall drop all 
the indices and simply write I{x) instead of 72 (aO and 
p(x,x,T) instead of p^(x,x,T). 



of the heavy particle, of mass M, with a light particle 
of the medium, of mass m. If M 3> m the time scale 
t s of a single-scattering process is short with respect to 
the typical time scale t of evolution of the heavy parti- 
cle. Thus, Born-Oppenheimer adiabatic approximation 
applies, and the dynamics of the center of mass of the 
heavy particle can be considered as frozen in the time t s . 
Accordingly, if 'J'(r) and x( r i) are ! respectively, the wave 
functions of the heavy particle and of the light particle 
before the collision, the wave function of the composite 
system after the collision is \l/(r)xr(r/), where Xr( r i) is 
the outgoing wave function of the light particle, scattered 
off at the point r = (x, y, z). As a consequence, the final 
state of the heavy particle is described by the reduced 
density matrix \l/(r) \l/(r')* (xr'IXr)- 

For arbitrary initial density matrix p(r,r'), and many 
independent individual scattering events, the variation in 
the time At of p(r,r') due to collisions is then 



Ap(r,r') ~ ~NAt[l-{ x Axr) J p(r,r'), 

where Af is the mean number of collisions per unit of 
time and (Xr'IXr) denotes the average with respect to a 
suitable ensemble of light particle wave functions. By 
taking into account also the rate of change of p due to 
the free dynamics, one arrives at the master equation of 
Joos and Zeh: 

^ = £ p + £ip, (18) 

where 
and 

(£ jP )(r,r')= -AAfl-(x7S7))p(r,r'). (19) 



A. Estimation of environmental coupling 



III. MARKOVIAN APPROXIMATION 

In order to evaluate 7(x), we need to determine 
p(r, r',i), the reduced density matrix at time t. In gen- 
eral, the evolution of p is highly non-Mar kovian, being 
the evolution induced by Eq. (|16fl via Eq. (|14|) . For an 
environment made of a gas at low pressure we may rely 
on the Markovian approximation provided by the model 
of Joos and Zeh Q . We shall now recall the essential in- 
gredients of this model and refer to literature for a thor- 
ough discussion 0, 0] . (For some basic steps towards 
a rigorous derivation see ; in this regard see also (2i| ; 
for a more general analysis of quantum Brownian motion 
see 112). 

This model aims at providing an autonomous evolution 
equation for an heavy particle moving in a gas of light 
particles under the approximation of negligible friction. 
To get a handle on the model, consider a single collision 



For a complete specification of the RHS of l|18fl . we 
need to evaluate the interaction operator H19|) related 
to the different processes inducing entanglement of Cgo 
with surrounding environment: scattering events (with 
thermal photons and air molecules) and photon emission. 
Such an evaluation of the interaction operator is rather 
standard, and can be found in the literature on the Joos 
and Zeh model — modulo some numerical values that we 
have corrected, and with the exception of our treatment 
of decoherence due to photon emission that is slightly 
different from what can be found in the literature (see, 
e.g., ii). 

1. Scattering with thermal photons 

In fullerene experiment, the environmental tempera- 
ture is 0£ « 300 K and thus the wavelength of thermal 
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TABLE I: Physical parameters 


of fullerene 


experiments |2|. 


Ivldbb (Jl ILlllelcllLi LvgO- 


M r. 


, 1 O Y 1 n-24 T7- 


Radius of C&q • 


D r~ 
IX ~ 


Q C v 1 A - 10 m 
O.O X lu 111 


_LLI1I1L)g1 dl Ul L (Jl 




~ yuu xv 


JZjIIV 11 (JllillL.il Ldl LLillJJL;! dLU.IL.. 




OUU IV 


Mean wavelength a of Cqq '. 


A w 


2.5 x 10~ 12 m 


Mean time of flight" 


T w 


6 x l(T 3 s 


Grating-screen distance: 


L = 


1.25 m 


Collimator aperture: 


a = 


1(T 5 m 


Effective slits width: 


d ps 


3.6 x 10" 8 m 


Grating period: 


D = 


l(T 7 m 



"Mean values are deduced by the measured velocity distribu- 
tion characterizing the fullerene beam outgoing from the oven (see 
Eq. 1611 below). 



photon is Ap/j = he/ (fcgQj) ps 4.8 x 10~ 5 m. As we shall 
see in Section llVl fEas. l|3fj|) - (|3"?|) and relative evaluation 
in Tab. Illlll . because of the incoherent preparation of the 
beam, we have that \x — x'\ < £q ps 10 _7 m. Under this 
condition, Eq. (|19|l assumes the form (see HE3) 



(£ I p)(x,x / ) 



»(scot; 



| x — x l \ 2 p(x, x') 



with 



^scat)_ g| 



C(9) 



2tt 



2.4 x 



io 2 urV 



( 2 °) 

where the fullerene molecule has been modeled as a di- 
electric sphere with the dielectric constant e r ps 4 pFj| 
and C(9) ~ 1, with £(z) representing the Riemann 
£ function. The previous relation has been obtained in 
the regime of Rayleigh scattering (since fullerene radius 
R ps 3.5 x 10~ 10 m is much smaller than X p h) and by 
using the Planck distribution for environmental photons 



2. Scattering with air molecules. 

Air molecules, with a mean mass m a i r ps4.8x 10 -26 Kg, 
at the temperature Qg « 300 K have a de Broglie wave- 
length X air = hj \Z2i:vn a i r kB^)s. ~ 10~ n m <C (-a (see 
Tab. GUI). Thus, assuming a Maxwell-Boltzmann distri- 
bution for air molecule velocity, it follows from the anal- 
ysis performed in [17J (Eq. (2.17)) that 



(£ip)(x,x') 



32tt 



k B @s m a 



for x 7^ x ' 



for x 



(21) 

where P{@s ) is the pressure at the temperature 6g and 
(Ttot is the total cross section for scattering events. 

In the case of fullerene experiment P ps 5 x 10~ 6 Pa 
and <T t ot ~ 9 x 10 _18 m 2 [54[. If F air (oo) is the constant 
value assumed by [C,ip){x,x') for x ^ x', from Eq. lj?T|) 



we get F a i r (oo) ps 32s -1 . In order to make a compar- 
ison between the different decoherence sources involved 
in diffraction experiments, we can introduce an effective 
localization factor A also for air scattering events. Given 
a pair of slits at the distance nD in a periodic grating of 
period D we have 



A air (n) 



■Fair(oo) 

(nD) 2 



(22) 



For adjacent slits (n = 1) the localization factor A a 
assumes its greatest value 



A air (n = 1) w 3.2 x 10 



15 . 



(23) 



Photon emission 



The model of Joos and Zeh can be extended to the 
description of photon emission processes. In fact, also in 
this case, the wave function of the composite system after 
a single emission event is, in general, \E , (r)x r (r;) I where 
*(r) is the initial wave function of fullerene and Xr( r i) 
is the outgoing wave function of the photon emitted in 
r. Since emission time scale, in analogy with scattering 
events, is much faster than characteristic time of fullerene 
free dynamics, the state |xr), in position representation 
and asymptotically in time, is well described by spherical 
waves 



{rt \Xr) oc 



,ik | rj— r | 



where k is the wave number of emitted photons. It follows 
that 



(Xr'IXr 



sin(fc |r — r '| 
fclr-r'l 



sinc(fc |r — r '|). 



According to Eq. (|19fl and assuming that fullerene diffrac- 
tion can be effectively treated as a one-dimensional prob- 
lem, the interaction operator for photon emission be- 
comes 



(dp)(x,x') = —M 



1 — sinc(fc | a; — x'\) 



p(x,x'). (24) 



In fullerene experiment there are essentially two chan- 
nels for photon emission: the black-body radiation and 
the disexcitation of internal vibrational energy levels. In 
particular, for black-body radiation at the fullerene tem- 
perature of Of ~ 900 K, we have a mean wavelength of 
emitted photon equal to 



A (Wj) 



1.6 x 10- a m. 



For decays of internal energy levels it was measured 
a peaked infrared spectrum with the shortest wave- 
length m 

XiZ b) « 7x HT 6 m. 
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In both cases it results (see Tab. 1111(1 A em ^> £q. This 
permits to simplify Eq. 12411 by the expansion of the sinc- 
function in powers of k\x — x'\. Keeping the first no- 
constant term, we straightly obtain 



(C I p)(x,x') = ~A em \x - x'\ 2 p(x,x'), 



with 



A. 



Nk 2 
6 



(in agreement with what obtained by Alicki |23j1. 

In the following, we shall calculate the mean value k 2 
for the two different channels of photon emission. 

a. Black-body radiation. The probability distribu- 
tion of the wave number k is given by the Planck law 



n(k) dk 



k R e f 



i 



k 2 dk 



2C(3) 



ft g jU 

,k B F K - I 



where e k 4.5 x 10 5 is the emissivity of fullerene at 
Of w 900 K [13 and C(3) w 1.2. Then Eq. becomes 



6 



3tt 2 C(5) We 



C(3) (AS') 2 



(25) 



with C(5) « 1.04. 

The number of emitted photons per unit of time 
can be estimated as JV ~ E^/(k B Q F ), where 
is the total energy emitted per unit of time and 
ksQ>F ~ 0.08 eV represents the single photon energy. 
By integrating the Planck distribution, one evaluates 
E^= eSa@% pa 16eV/s, where S = knR 2 is the total 
surface of fullerene macromolecules (R pa 3.5 x 10~ 10 m) 
and a is the Stefan-Boltzmann constant. It results 
TV pa 200 coll/sec, and thus 



energy levels for fullerene, we directly refer to the results 
of experimental spectroscopy. Since infrared spectrum 
shows a peaked structure, we can write 



A 



(vib) 

P h 



^<®!, 5 xl0W, (27) 



G 



(i 



where k* pa 9 x 10 5 m is the wave number related to 
the most energetic spectral line (see Fig. 4 in [2^) and 
TV pa 400coll/s [J. 

A direct comparison between evaluation l)2()[l. 123(1 . 
(|26|l and (|27() . reported in Tab. UTI shows that the main 
decoherence processes are scattering with air molecules 
(especially for adjacent slits, cf. Eq. iJUJl), followed by 
photon emission due to decay of internal vibrational en- 
ergy levels. 

TABLE II: Sources of decoherence in the conditions of 
fullerene experiments Q, Q. 

Decohering event A (m 2 s _1 ) 



Scattering 

with thermal photons: 

with air molecules: 

Photon emission 

black-body radiation: 

decay of excited states: 
Global effect: 



A (scat) ^ 2A x 1Q 2 

Kir < 3.2 x 10 15 . 



Pd 2.5 x 10 9 , 
A^ b) < 5 x 10 13 . 

A < 3.3 x 10 15 . 



B. The effective master equation 

According to the foregoing analysis, the Joos-Zeh equa- 
tion l(T5jl effectively reduces to 



A 



(66) 
ph 



2.5 x 10 9 m 



(26) 



b. Decay of internal vibrational energy levels. Since 
we lack of a model able to describe decays of internal 



ih 



dp 
dt 



[H 0l p] - iA[x, [x,p]] 



or, more explicitly, 



ih dpix^_ = t?_(£ 5 _0^ p(aj x , it) _ lA {x _ xTp{Xj ^ t) (28) 

with 

A = A air + A ph , 

where A a i r is given by Eq. ill' I'll and the three terms in A p h = Ap™ t ' ) + + A^ are given respectively by 

Eqs. (EOJ, and (S7jl. 

Eq. (|28[1 is a well-known equation and its solutions are readily obtained (see, for example, App. 2 of Joos in |l8j ) 

p(x,x',t) =11 dx dx' K(x,x',t;xo,x' o ,0) po{x , x' ) , (29) 
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where 

(x - x') 2 + (x - x' ) 2 + (x - x 1 ) (x - x' Q ) I . 

(30) 

Notice that the first exponential describes the free dynamics, while the second takes into account the interaction with 
the external environment. , 



K(x,x' ,t;x ,x' ,0) = 




IV. PREPARATION OF THE INITIAL STATE 

In order to determine p(x, x', t) and thereby evaluating 
the intensity on the screen given by Eq. 117(1 , we still need 
to specify the initial density matrix po(x,x), taking the 
initial time t = at the moment of the splitting produced 
by the diffractive grating. 

Because of thermal production and in spite of the fol- 
lowing collimation, each fullerene wave function has a 
(mean) transversal wave number k x (ideal collimation 
would correspond to k x = 0). Thus, after the splitting, 
the macromolecule wave function is of the form 



(31) 



where tp s represents the s-th of the N slit-shaped wave 
packets outgoing from the grating. The beam is an inco- 
herent mixture of such wave functions with wave number 
k x randomly distributed according to a probability dis- 
tribution p (k x ). This distribution depends on the geom- 
etry characterizing the collimation setup, which reduces 
the wide thermally produced spread on x direction. The 
density matrix of the beam at t — is then 



p (x,x 



dk x p(k x ) ip (x;k x ) ip (x';k x ) 



Letting 



po(x,x') = y^y s (x) ip s '(x'Y 



(32) 



(33) 



we obtain 

p (x,x') = p (x, x') J dk xP (k x )e- lkAx '- x) 

= V2Trpo(x,x')p(x' - x) , 

where p is the Fourier transform of p and po(x, x') has the 
meaning of the density matrix in the ideal case of per- 
fect collimation. A typical diffraction setup consists in a 
periodic grating of period D, which we consider placed 
symmetrically with respect to the optical axis as illus- 
trated in Fig. |5J). In other words, we consider 



Po 



(x,x') = J2 ( 



D 



lfi[X + S— ] lfi[X + s — 



,D 



(34) 



where s,s' = ±1, ±3, ±(N — 1) (for symmetry with 
respect to the optical axis, N is considered to be even). 
The size Ax of the support of p is simply fixed by 



The general structure of Eq. Li-it (for a treatment of 
which we remind also to the section 9.1 of Joos in [Is|^ 
appears for any choice of the density matrix po(x, x') and 
in every case in which a particle is subjected to an uncon- 
trollable source of random "kicks" which produces instan- 
taneous shifts in momentum, as in Eq. I|31[) . Moreover, 
in case of random kicks with a mean momentum transfer 
position-dependent (for example, in case of van der Waals 
interaction between crossing particles and atoms of the 
grating) , the effect on the initial state consists in an effec- 
tive reduction of the aperture width . A similar effect 
in molecular diffraction has been already investigated in 
the framework of classical optics [2|J . 

Now, in order to simplify the analysis, we shall adopt 
the convenient and physically reasonable assumption of 
a Gaussian probability distribution 



P(k x ) 



1 



exp 



2 < 



so that 

p (x,x') = p (x,x') exp 
where we have defined 



(x - x' f 



211 



(36) 



(37) 



Ax - ND . 



(35) 



This quantity, which will play a relevant role in the 
following analysis, will be called the coherence length (at 
time t = 0). We note that for £q < Ax there is a bound 
on the length on which the macromolecules can be co- 
herent, expressed by £q itself. In particular, only beams 
characterized by an initial coherence length £q > D may 
produce a coherent superposition of wave packets, and 
thus interference fringes, on the detection screen. On the 
other hand, for £q ^> Ax, the damping shown in Eq. I|3t)|) 
does not take effect, and the preparation of the initial 
state results to be coherent on the whole support Ax. 

Now we consider explicitly a typical diffraction exper- 
iment with macromolecules (see Fig. |2J), in which the 
collimation apparatus consists in two identical slits with 
aperture a, at a distance I S> a. The greatest drift ve- 
locity along the x direction results to be |wj;| m ax = v6, 
where v is the macromolecule classical velocity along the 
optical axis and 8 = a/l is the angle under which a point 
situated in the aperture of the first collimator sees the 
aperture of the second one (since I 3> a, then the angle 
9 can be considered the same for every point of the first 
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Collimation 



Diffraction 




+ D 



FIG. 2: Collimation setup and diffraction grating. For clarity 
the diagram is not in scale. 



collimator). Thus we can put 3<7fc x = MvO/h and so we 
obtain 



Ml 



[Mi 



V 3h 



(38) 



An evaluation of the initial coherence length £q for 
fullerene experiments is reported in Tab. IIIII 



V. THE INTERFERENCE PATTERN 

Consider now the initial density matrix po given by 
Eq. (O- Define 



1 



2£(t) 



At 1 
~3~ 



to 



2At p2 
3 Z 



(39) 



(note that £(0) = £ ). Then Eq. JDJ becomes 



p(x,x',t) = 



(x-x ) 2 -(x'-x' )' 



(x - x' Q f At 
2£(tf 3 



(x-x') 2 + (x-x')(x -x' ) 



whence, from relation l|17|) . the intensity on the screen is given by 



h(x) = p(x,x,T) 



m f f t r %m 

Axq dx exp 



2nhT 



2hT 



x a - x o + 2.x(a;Q - x ) 
I 



Po(a;o,^o), 



2P 



pv(x ,x' a ), (40) 



where xo and Xq run along the slits crossed by the initial 
wave function whose support is Ax and £ = £(T) is the 
coherence length at the time of flight t = T. 

As already argued for £q, Ml < Ax the exponential 
exp[— (x ~ x' ) 2 /(2£ 2 )] reduces from Ax to £ the length 
scale on which the initial state is coherent. This scale 
is fixed from both geometry of the experimental setup, 
i.e., the collimation apparatus and the distance L be- 
tween grating and screen, and the physical conditions 
under which the interferometry takes place, i.e., the mo- 
mentum of the macromolecule and the effect of the en- 
vironment (see Eqs. I|39|l and JSJ)). On the other hand, 
the above exponential does not give any relevant contri- 
bution if £ > Ax, and, from Eqs. and it 
follows that this occurs when there are both good colli- 
mation (9 w 0) and negligible coupling with surrounding 
environment (A as 0). In this case we fall back to the 
treatment of section III Al 

Note that interference fringes appear on the detection 
screen if Ax > D and £ > D, i.e., if the molecules are 
coherent at least on two contiguous slits. The numerical 
estimates of £q and £ in the condition of fullerene exper- 
iments are shown in Tab. IIIII In particular, note that 
£ D and thus interference is mainly due to adjacent 



slits. Moreover, a comparison between the values of £o 
and £ shows that the main mechanism which yields a loss 
of coherence is the angular divergence of the beam 55] . 



TABLE III: Comparison between losses of coherence in 
fullerene experiments Q, Q. 

1.3 x l(T 7 m 



Initial coherence length (t = 0): 
Coherence length at t = T: 



£o 



D = 10" 



A. Fraunhofer approximation for the intensity 

We shall now proceed to an approximate evaluation of 
I(x), relying on conditions that are reasonable in com- 
mon interferometry experiments performed in far-field 
approximation (see App.^for a more refined evaluation 
in the case of a pair of Gaussian shaped slits). 

Note that 



exp 



iM(a 



2hT 



10 



when 



M(xl - x' a 2 ) _ M(x + x' )(x - x' ) 



2hT 



2HT 



< 1 



and this condition is clearly satisfied in the Fraunhofer 
regime (|13f) . Nevertheless, in the presence of a coher- 
ence length t < Ax, we have relevant contributions in 
integration (|40l) just for (xq — x' ) < I. Thus, under the 
condition 



MAxl 

hT 

l2(x) gets approximated by 
M 



< 1 



h(x) 



2nhT 



dxo dx' exp 



iMx 



{x' -Xq) 



(41) 



pi(xo,x' ), 
(42) 



where 



Pi(xo,x' ) 



exp 



(x - x' ) 2 



2P 



p (x ,x' ) . 



Introducing the Fourier transform of p\ 
we have 

T , . M ^ T r — Mx 2irx , A _ S 

h{x) = —pi{k,-k), where k = — = — . (43) 

This result is very analogous to Eq. of section III Al 
with (|41|) replacing (|1U|) whenever t < Ax. Also in this 
case l|41|) should be rewritten in terms of the physical 
variables under control (cf. (|12|l 'l. namely as 

t_ 2i_ 

L ^ ~Ax ' 



This notwithstanding, there are some basic difference 
that should be underlined: First, p\ is not the initial 
state, but it is an effective state that takes into account 
incoherence due to preparation and to the time evolution. 
In fact I depends on the physical and geometrical vari- 
ables of the experiment in the phase of preparation and in 
its future development and it is progressively reduced by 
increasing the time of flight. Second, unlike what typ- 
ically happens in the framework of classical optics and 
the theory of scattering, it is no more useful to evaluate 
pi asymptotically in time, since environmental-induced 
decoherence completely destroys interference fringes at 
times too large. 



In fullerene experiments, 



10 



-7 



< Aa 



(for an estimate of Ax, see Section FVlip . In this case the 
left-hand side (LHS) of l|41|) is not at all negligible with 
respect to unity. Anyhow, a more precise inspection of 
integration H40|) , with po given by l|34(l , shows that condi- 
tion (|41|) is a too strong demand and that approximation 
(|42|l can be reasonably applied. So doing, the error made 
is not completely negligible only for the pair of adjacent 
slits farthest with respect to the optical axis. This error, 
however, affects negligibly the sum involving the contri- 
butions of all the slits. 



Let us now compute pi(k, —k) for po expressed by <|34[1 . 
i.e., for the split of macromolecules on a periodic grating 
of period D. First we make the change of variables £ = 



x Q + sD/2 and f 



s'D/2, which leads to 



Pi 



3 ik(s — s') 



D [( s -s')P/2] 2 



2 g 2£ 2 



<p(o v>(t'r , (44) 



and second we perform the Taylor series expansion of the real exponential in the previous integral in the variable 
(C - It and about the point £' = £ 



exp 



2^ 2 



2P- 



(45) 



The solution of Eq. (|44|l is particularly handy whenever 
the effects due to incoherence are negligible on a length 
scale of the order of the slit width d or, in other words, 
whenever the strength of incoherence does not spatially 
resolve the single slit. This is assured by a slit width 



much less than the coherence length 

d/l<£l. (46) 

Under condition (|4(j|l the LHS of Eq. I|45|) is well approx- 
imated by 1. In fact, since £ and £' run within the slit 
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width d and (s' — s)D/2 < I due to the damping ex- Thus the Fourier transform Q44[) becomes 
ponential exp j- [(s - s') £) / 2 ] 2 /( 2 ^ 2 )} in E q- dSJ, we 



have 



2£ 2 



pi(k,-k) 



2tt ^-f 



[( a -s')P/2] 2 



r 



(Although the rough condition 14(i|) is not directly satis- 
fied in fullerene experiments, the zero-order approxima- 
tion of (|45|l can be reasonably applied in integration (|44H ; 
for an evaluation of the error introduced the reader can 
see App.|Bj). 

In the light of these considerations, the intensity pat- 
tern is well approximated by 



h{x) 



M 
KT 



where is the Fourier transform of ip. 



Note that the sum of the terms with s = s' simply 
jives N, while the sum of the terms with s s' gives 



E [(s'-s)D/2] 2 
e i? 2 



cos 



k{s' - s) 



By adding these two contributions and with 
k = 2irx/(\L) (see Eq. we arrive at the suggestive 

"Fraunhofcr-like" expression 



h(x) = I{x) 



2tt 



^(2-kx 



N- 



N-l 



(N-n)e ~^~cos 



71=1 



1 XL 



(47) 



r 



(being understood that for N = 1 the sum is zero) , where 
n=(s' - s)/2 0. 

Equation (|47|l shows that, whereas all the N wave 
packets outgoing from the grating contribute to the in- 
tensity revealed on the screen, the pairs of slits which 
concretely contribute to interference oscillations are dis- 
tant at most of the order of £, due to the damping ex- 
ponential in the sum. It follows that, for a finite £, the 
interference pattern shows "distortions" in fringe struc- 
ture due to partially random preparation and decoher- 
ence, but, being incoherent effects typically negligible on 
single-slit space scale, fringe pattern is just modulated by 
the single-slit diffraction profile, |(J5| 2 , according to clas- 
sical optics Fraunhofcr diffraction. 

As already sketched before, it should be observed that 
the intensity on the screen may show interference fringes 
only if the coherence length is at least as long as the 
grating period, i.e., 

£(t) > D (48) 
(note that this inequality should be satisfied at least at 



the initial time t = 0). For positive times, recalling 13911 . 
we obtain 



t < T P . = 



2AD 2 



which provides an upper bound for the time of flight, i.e., 
an evaluation for the effective coherence time r c . Note 
that for the fullerene experiment, according to Tab.[I]and 
Tab. Hn it results 



r r = 4 x 10 _2 s . 



which is indeed several times the value of the time of 
flight in this experiment |57j . 

The effective coherence time t c is clearly an upper 
bound for the time of flight T, since interference fringes 
are detectable only within t c . This shows that the "geo- 
metrical optics" limit in presence of decoherence requires 
more care than in the free case. In particular, it can not 
be based on the standard time-independent methods and 
the ll t — > oo" limit. 
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B. Extension to a generic angular divergence of the beam 

This section is devoted to generalize Eq. 1(47(1 for a generic transversal wave number probability distribution p (k x ). 
Introducing Eqs. ((30(1 and ((33|) in Eq. 1)29(1 . the long-time asymptotic behavior of the intensity becomes 



h(x) 



M 



By the same variable change which leads to Eq. ((44(1 and making explicit po for a grating of period D (cf. Eq. 1(34(1 ). 
we obtain 



M 



27r^T 



e -^f [(s'-s)D/2] 2 



(49) 



As discussed in the preceding section, note that |£' — £| < d and that |(s' — s)D/2| max ~ (AT) x / 2 due to the damping 
term e~^r^ s ~ a ) D / 2 ] in Eq. ((49(1 . Thus, in case of decoherence negligible on the single-slit length scale, i.e., for 
d(AT) 1/2 < 1, and for a slowly varying function p, such that p[± d+ (s - s')D/2] ~ p [(s' — s)D/2] Vs,s', Eq. 
becomes 



I^tt) 3 / 2 



27ra; 
AT 



iV + 2^(A-n)e~ 



p(nD) 



p(0) 



2-kDx 



(50) 



where n — (s' — s)/2 and fc = 2t:x/(XL). Note that the assumption of slow variation of p is directly assured by a 
sufficient sharpness of the wave number distribution p(k x ), i.e., by Ak x -C d . 



r 



VI. QUANTUM INTERFEROMETRY AND 
CLASSICAL DIFFRACTION THEORY 



A. Comparison with geometrical optics 

In case of complete coherence, i.e., for p{k x ) — S(k x ) 
and A = 0, Eq. 1(50(1 reduces to the well-known Fraun- 
hofer relation for optical diffractive patterns 



/(*) 



2tt 
XL 

2?r 
XL 



N-l 



Kf ) E 



2tvDx 
1 XL 



2ixx 
AL 



n=0 
2 



sin(7rA^L>x/AL) 
sin(7rDx/AL) 



Moreover, for just two slits (N = 2), i.e., for Young 
double-slit interference, Eq. 1(5011 becomes 



I(x) 



47T 

XL 



P(0) 



2nx 
~XL 



p(jlD) AM XL 

1 H — e 



p(0) 



6"" cos 



(2-kDx 
1 AL 



This expression is very similar to that used in classical 
optics to describe interference patterns due to partially 
coherent electromagnetic fields ^ij. I n particular, note 
that the damping term for quantum interference oscilla- 



tions 



3M 



p(nD) KM XL 

— i p 6jtR j 

p(Q) 



(51) 



is the quantum-mechanical counterpart of the fringe vis- 
ibility Vco °f classical optics 



Vqm 



Vco- 



(52) 



Both for quantum and classical interferometry, the visi- 
bility V is a measure of the distinctness of the fringes and 
is defined by 



V = 



-MTiax -Mr 



(53) 



The intensities / max and J m i n are, respectively, the maxi- 
mum and the minimum revealed on the detection screen 
in the immediate neighborhood of the optical axis. 

Now it is useful to recall an important result from the 
classical theory of partial coherence. The pattern visibil- 
ity Vco of a quasimonochromatic field, equally split by 
a pair of slits, coincides with the modulus of the spectral 
degree of coherence /i(A) [I^Eul]. which characterizes the 
field correlation in the space-frequency domain 



V, 



CO 



KA)|. 



(54) 



From Eqs. ((53(1 and 1(54(1 it follows that the degree of spec- 
tral coherence is upper bounded by unit, value assumed 
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in condition of complete coherence (e.g., in case of laser 
radiation diffraction) . 

According to the correspondence l|52|l . the results of 
the classical theory of partial coherence extends to quan- 
tum systems, mutata mutandis. For instance, in Sec- 
tion ^^H] we shall show some interesting analogies con- 
cerning with temporal and spatial coherence of beams, 
while in the following we underline the differences exist- 
ing between classical optics and quantum mechanics. 

First of all, the degree of coherence of quantum parti- 
cles depends both on the collimation of the macromolec- 
ular beam and on the strength of interaction with the 
surrounding environment during the time of flight (cf. 
Eq. f5TJl ). In optics, instead, the degree of coherence 
of quasimonochromatic fields is only due to source de- 
tails. More particularly, the corruption of visibility of 
interfering fields increases with the spatial extension of 
the source, composed by a statistical ensemble of many 
independent elementary radiators. 

Moreover, in classical optics the explicit form of the de- 
gree of coherence /z(A) depends on the geometrical shape 
of the source, while the damping term Vqm depends both 
on the features of the evolution kernel (|30fl , characteriz- 
ing the decoherence model, and on the geometrical details 
of the collimation apparatus (cf. Eq. (JS^)- 

B. Fresnel regime and Talbot interferometry 

We would like now to comment on interferometry in 
the near-field zone [33, 113 , which has been re- 

cently realized by means of Cm beams [3(| [13 ■ Such 
experiments show that, at distances from the diffraction 
grating multiple of the length 2D 2 /X, images of the grat- 
ing itself are reconstructed (see also the optical Talbot 
effect HHil). 

Thus, by shifting another identical grating, placed be- 
hind the previous one at a distance 2D 2 /X, the integrated 
signal outgoing from the gratings periodically changes 
from its minimum (half period displacement of the two 
gratings) to its maximum (complete alignment). 

If the influence of the environment is negligible, a treat- 
ment of this effect in the spirit of Section III Al can 
be performed. In fact, given the correspondence be- 



tween Helmholtz and stationary Schrodinger equation, 
one can directly exploit the standard optical techniques, 
such as Fresnel-Kirchhoff diffraction integrals in Fresnel 
zone, with suitable boundary conditions — "transmission 
functions" — at the gratings. Indeed, this is what it has 
been done (see, e.g., [lEEl) by means of the so-called 
"paraxial approximation," assuming both gratings dis- 
tances large with respect to the grating period and an 
infinite number of slits. 

In experiments with large molecules [37]], it has been 
observed that the visibility of the signal is progressively 
reduced by increasing the pressure of environmental gas, 
a clear sign of environmental quantum decoherence. A 
quantitative explanation of this effect — using the model 
of Joos and Zeh in order to suitably modify the classical 
Fresnel-Kirchhoff description recalled above — has been 
already provided in |37| . A more self-contained and thor- 
ough analysis, based on Eq. 1|17[1 . will be presented else- 
where [23. Here we shall provide just a sketchy outline, 
referrin g to a theoretical treatment already present in lit- 
erature [4j . In this last work, the form of the propagator 
describing the free evolution of a quantum wave is shown, 
split by a diffraction grating with a formally infinite num- 
ber of slits, i.e., with an associated initial density matrix 
of the form 

Po(x,x') = ^2 f( x + 3 D ) p( x ' + j' D T ( 55 ) 

(note that, with respect to initial state expressed by l|33|) 
and (|34|) . here it is assumed perfect collimation and an 
infinite number of slits). Let K^ rcc ^ be the so called 
Talbot propagator, including the sum over j and f of l|55|) 
and, accordingly, providing the intensity pattern I(x) in 
term of the single wave packet f(x): 

I(x) cx J dxo dx' K^ ICC \x, x,t; xo, x' ,0) ip(xo) ip(x' )* . 

l-l (56) 
It results (see e.g. |4l|) that, at a distance Lt from 

the grating equal to 2D 2 /X, or multiple of it (and con- 
sequently at times tr multiples of Lt/v = 2MD 2 /h), 
K}p Tee ^ reduces to 



K^ ICC \x,x,t T ;x o ,x' Ol 0) = exp 

j=—oa 
+00 

- J2 6(x-x +jD) J2 S(x-x' 0+3 'D). 

j = ~co j'=— oo 

Clearly, from this relation and Eq. H5(jf) it is immediate to verify that the final state is an exact reconstruction of the 
initial one (|55|l . 

By following the same steps leading to Eq. (|57|l , but now using the propagator l|3U|) which embodies the incoherence 



2nij- 



•i'O 



D 



+ OC 



exp I — 2irij' 



D 



(57) 
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effects, the Talbot propagator becomes 



K^ nv \x,x,t T ;x 0l x' ,Q) = —2 ^2 exp(-2irij— — ^ ) exp <J - 2j 2 



exp I - 27rr/' 



J'=— oo 



exp ^ - (2j" + 4j 3 >) 



D 



where i(tr) is the coherence length computed at the Talbot time tr- With respect to Eq. (|57|l . here some additional 
exponential are present which obstructs the reconstruction of the initial wave function. So, if we put a second grating 
at distances multiples of 2D 2 /A, even for a perfect alignment between the two gratings the wave function is partially 
stopped and thus the intensity detected will be lower than in the free case. Similarly, for a displacement between 
the gratings of an half period, a portion of signal, even little, may be detectable further on. In such a scenario, a 
decrease of the coherence length £, i.e., a growth of the incoherence of the beam, leads to a progressive reduction of 
the visibility of the total intensity on the screen, in agreement with the behaviour of the experimental data reported 
in [371 . (An improvement of this analysis should presumably take into account: 1) the three free-standing gratings for 
the Talbot-Lau interferometer effectively used in experiments; 2) a proper description of the van der Waals interaction 
with the grating). 



C. Near- field interferometry and randomness of 
arrival times 



Near-field interferometry, such as Talbot-Lau interfer- 
ometry, should allow us to probe quantum effects also due 
to the motion along the longitudinal direction, which so 
far has been treated as classical. Such a treatment has 
been of course completely motivated by the experimental 
conditions considered so far for which both (J5J and Q 
are satisfied with a high degree of approximation. But 
suppose that position spread in the longitudinal direc- 
tion Ay is not completely negligible with respect to L. 
Then the arrival times would have statistical fluctuations 
of order 

AT - Ay /v. 

For the fullcrenc experiments in the Fraunhofer region 
such fluctuations are not appreciable: in this case, 



AT ~ Ay/v 
we have 



5 x 1CT 9 s, and since T w 6 x 1CT 3 



s, 



AT < T. 



(58) 



Near-field interferometry, possibly with light particles, 
should be able to test the measured intensity when (|5H|> is 
violated. A first prediction is immediately suggested by 
©: the measured intensity is obtained by the intensity 
distribution consider so far, namely p(x, x,T), by convo- 
lution with \4>o\ 2 - Thus, randomness of the arrival times 
appears as an independent noise on the standard inter- 
ference profile, which reduces the fringe visibility as it 
were an additional source of "decoherence" . This effect 
could be confused with a sort of intrinsic decoherence 
(in this regard see a recent proposal concerning atomic 
diffraction by standing light wave |42^1. 

More generally, one may analyze the predictions of @ 
in the mesoscopic regime []] . Let us consider a monochro- 
matic beam, devoid of angular divergence, composed by 



free light particles of mass to. Let the beam be diffracted 
by two slits of width d and distance D. It is convenient 
to describe the split wave function ^(r) by means of two- 
bidimensional Gaussian wave packets, whose barycen- 
ters move parallel to the optical axis y with the velocity 
v = hky/m: 



*(r) = i> (x)My) 



= c 



iD/2) / 





L s=± 



■ + ik v y 



(59) 



The constant C ensures normalization. The transversal 
standard deviation <r x of each wave packet is related to 
the slit width d (a typical assumption is &a x — d), while 
the longitudinal one is related to the extension of the 
initial support of ^(r). 

After the splitting, we can consider that Vf'(r) evolves 
according to the free Schrodinger equation J2J) until the 
particle is detected on the detection screen placed at a 
distance L. Given an arbitrary distance, not necessar- 
ily large compared with the size of the support of \&(r), 
and for a free evolution, the intensity detected at y = L 
can be calculated by means of Eq. Q, with ip(x,t) and 
4>{x, t) obtained by the free evolution of the initial state 
(|59|l . The result of such a numerical simulation for ultra- 
cold neutrons is shown in Fig. [3J At distances as short 
as to be comparable with the longitudinal spread of neu- 
tron wave packets, the statistical fluctuations on arrival 
times produce an appreciable reduction of the fringe vis- 
ibility just as it would happen in case of an incoherent 
preparation of the beam and/or in case of environmental- 
induced decoherence. The only difference consisting in 
the distance dependence of the different processes: this 
kind of "decoherence" reduces by increasing the distance, 
environmental-induced decoherence increases, while the 
effects due to incoherent preparation are independent. 

In particular, notice that the reduced visibility is not 
due, not even partially, to an incomplete wave-packet 
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V = 0.70 



V = 0.96 




Position (mm) 



FIG. 3: Double-slit interference patterns for different values 
of the distance between the grating and the detection screen. 
Patterns are obtained for ultra-cold neutrons (v m 1 m/s) 
diffracted by a pair of slits of width d — 5 x 10 m and sep- 
arated by a distance D = 10 -5 m. Neutron initial wave func- 
tion is described by means of bidimensional Gaussian wave 
packets with a longitudinal standard deviation a y — 10~ 3 m 
and a transversal one o x = d/6. 



distribution of a strictly monochromatic beam. Some ad- 
justments have to be carried out in order to include in 
our treatment the effects on the diffraction pattern due 
both to the velocity distribution characterizing the beam 
macromolecules and to the distortions unavoidably in- 
troduced during the measurement process |45| . Each of 
these corrections have to be implemented on the intensity 
level, since they consist in incoherent contributions. 

The total intensity is obtained by the sum of the 
monochromatic components of the beam 



I(x) 



d\f(X)I(x,X), 



(60) 



where I(x, A) is given by Eq. l|T7Jl and its dependence on A 
is shown by Eq. ©. The wavelength distribution /(A)dA 
is directly obtained by the supersonic velocity distribution 
f(v)dv characterizing the macromolecule ensemble |2j 



f(v) dv ex v 3 exp [ — (v — vq) 2 /v 2 ~\ dv . 



(61) 



It describes beams in transition between effusive and slow 
jet sources 46f. The parameters v$ and v depend both 
on the temperature of the beam and on the physical fea- 
tures of the given molecule and they are deduced by a 
best fit over experimental measurement of the velocity 
distribution (see Fig. 2 in |2(). 



superposition, being this ensured, also for the shortest 
distance shown in Fig. [21 by the large wavelength of ultra- 
cold neutrons. Thus, the loss of fringe contrast has to 
be ascribed only to arrival time fluctuations of the same 
order of the classical time of flight. 

For simulations of Fig. [31 we used a longitudinal der- 
ealization at the double-slit given by a y = 10~ 3 m. This 
assumption can be relaxed to shorter values still detect- 
ing results analogous to those shown by Fig. [3J but at 
shorter distances. In this case the spatial resolution re- 
quired for efficient detection becomes higher. Conversely, 
for a larger longitudinal derealization, intrinsic decoher- 
ence effects can be readily detected at larger distances by 
means of less refined detectors. 

In conclusion, let us underline that, since very slow 
neutrons (v « 1 m/s) 1 13l characterized by a wide 
transversal support (Ax > 10~ 4 m) ^3|, 0] have been 
used in intcrfcromctry, we think that an experimental 
test of the predicted behavior shown in Fig. [3] might be 
indeed in the reach of present technology (not necessarily 
for the concrete situation we have simulated, which was 
mainly for illustrative purposes). 



VII. NUMERICAL CALCULATIONS 
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-150 -75 
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FIG. 4: Theoretical fit and experimental data for Ceo fn- 
terferometry 0- Theoretical fit was performed consider- 
ing N = 10 equal, rectangular shaped, slits of an effective 
width d — 36 nm. The spatial resolution of the detector is 
2x = 8 /im. The theoretical curve and experimental data are 
normalized to the value of the central maximum. 



In this section we shall use Eq. (|47|l to fit the ex- 
perimental data reported in In this regard, note 
that Eq. (|47|) describes an ideal situation where an in- 
finitely accurate detector measures the spatial intensity 



The finite spatial resolution can be taken into account 
by an integration over the size of the elementary detector, 
say 2x, weighed on its spatial response function D(x). In 
particular, in the case of a flat response function, D(x) 
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becomes the characteristic function defined on the in- 
terval [ — x , +x\ and the effective detected intensity is 
expressed by a moving average 

I e s(x) = [d(I(x + ()D(0 = i= [7(1(0, 



whence, from Eq. 1)60(1 . 




150 



Position (/im) 

FIG. 5: Interference patterns due both to a macromolecu- 
lar beam characterized by a velocity distribution given by 
Eq. IbH as in |j| (full line) and to a strictly monochromatic 
beam corresponding to the mean velocity (dashed line). The 
curves are normalized to the value of the central maximum. 

For the free gap d of the rectangular slits, it was 
used, according to [2|j, the effective width estimated in 
Refs. pi as previously discussed in Section llVl 58] . 
Moreover, a constant background has been subtracted 
from experimental data 0, • 

Finally, referring to Eq. 147|) , note that the number N 
is basically a priori unknown, being the effective num- 
ber of slits which concretely contribute to diffraction. In 
fact, only the knowledge of the fullerene wave function 
on the grating would allow us to fix the effective number 
of slits in which the initial wave function is split, since 
the width of collimation only gives information about 
the maximum number of wave packets which could con- 
tribute to the interference pattern. Nevertheless, N can 
be easily inferred a posteriori as a free parameters of the 
fit with experimental data. So doing, we find N = 10. 

Our theoretical evaluation for the interference pattern 
is shown in Fig.QJ together with the experimental points 
published in Ref. 0, page 2819, Fig. 5. 




FIG. 6: Diffraction of molecules with masses multiple of the 
Ceo mass M. For larger masses the quantum behavior be- 
comes progressively negligible, approaching to the classical 
limit. Every curve has been normalized to the value of the 
central maximum obtained for the mass M. 




FIG. 7: Interference patterns for different values of the dis- 
tance between the grating and the detection screen, in unit of 
the value L reported in Increasing the distance the pat- 
tern spreads in position with respect to the optical axis. The 
visibility is reduced by the increased number of decohering 
events. Every curve has been normalized to the value of the 
central maximum obtained for the distance L. 



VIII. CONCLUSIONS AND PERSPECTIVES 

One of the main goals of our work has been that of 
deducing, within the framework of Joos and Zeh model, 
a geometric optics limit of quantum mechanics in the 
presence of decoherence. 

Our theoretical analysis confirms the negligibility of 
environmental disturbances in recent experiments of 
macromolecular interferometry 0, Q with respect to the 
loss of coherence due to beam production. 

Moreover, our analysis reproduces classical results, as 
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FIG. 8: The two figures show interference patterns for increasing pressures at fixed temperature, in unit of the pressure 
P = 5 x 10 _6 Pa used in 0- Left figure corresponds to the experimental setup reported in 0: slit effective width d — 3.6 x 10 -8 m 
and collimator aperture a = 10 -5 m. Right figure has been obtained for narrower slits (d = 2 x 1CF 8 m), an improved collimation 
(a — 5x 10 -6 m) and a velocity selected 6*60 beam (Av/v = 10%, where v is the mean velocity). Every curve has been normalized 
to the value of the central maximum at the pressure P. 



well as Fraunhofer relation for optical diffractive pat- 
terns, and provides the quantum-mechanical analog for 
interferometry with partially coherent sources of radia- 
tion. For the latter subject, referring to the theory of 
partial coherence, existing analogies can be pointed out in 
a deeper way by exploiting numerical simulations. Fig. [SI 
clearly shows that, by selecting the particle beam veloc- 
ity, the visibility of the corresponding interference pat- 
tern does not tend to unity, being upper bounded by the 
damping term Vqm (cf. Eqs. lj5Tl) - (|55)l ). On the other 
hand, the effect of the velocity selection makes more in- 
terference fringes visible at the border of the interference 
pattern. 

The same behavior is obtained in the framework of 
classical optics, studying interference patterns due to 
quasimonochromatic fields [4?1 |4S| . Classical partial co- 
herence theory, supported by recent experiments 
states that, by filtering, the pattern visibility at most 
approaches the value of the modulus of the degree of 
spectral coherence /u(A), which depends on the spatial co- 
herence of the source, and that more fringes becomes 
visible, since the temporal coherence is improved. 

Moreover, our approach provides a useful theoretical 
framework for analyzing present (and possibly new) in- 
terference experiments. For instance, we have studied the 
mass dependence of the interference pattern due to the 
diffraction of heavy particles. Fig. El shows the simula- 
tions corresponding to beams of macromolecules heav- 
ier than Ceo, but characterized by the same physical 
features. Note that for larger masses the quantum be- 
havior becomes progressively negligible, approaching to 
the classical limit. There are, indeed, other experimen- 
tal investigations which support the previous expected 
behavior [53. They show that visibility of C70 diffracted 
beams is slightly reduced than that obtained for Cgo ones 
in the same conditions. 



In studying the effects of decoherence, it is particu- 
larly interesting to analyze the case of a diffracted quan- 
tum particle which experiences a lot of scattering events 
before reaching the detection screen. This situation is re- 
alized for larger grating-screen distances (see Fig. Q) , or 
for more frequent scattering processes due, for example, 
to increasing value of pressure (see Fig.[SJ. In both cases 
the fringe visibility, and thus the wavelike behavior of the 
molecule, is progressively corrupted. 

By improving the collimation of velocity selected Cgo 
beams, the interference pattern shows a richer structure 
of fringes and thus a more evident quantum behavior. 
Moreover, interference oscillations appear also in less 
restrictive environmental conditions, provided that the 
signal-to-noise ratio is such to allow experimental detec- 
tion. In fact, as shown by the right plot in Fig. side 
maxima, far from the optical axis, which are not detected 
at pressure P in experiments P, Q , turn out to be clearly 
visible even for pressures ten times larger than P. This 
might be relevant in devising new interferometry experi- 
ments directed to the study of the quantum behavior of 
macroscopic objects and also to test quantitatively the 
effects on a quantum subsystem due to external noise. 
Our suggestion has been partially realized in very recent 
experiments |5lj . which show results in agreement with 
the prediction of Fig. [H] 
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APPENDIX A: YOUNG INTERFERENCE PATTERN WITH GAUSSIAN SLITS 

In the following we shall develop the exact solution of (|40|) in the case of slits with a Gaussian shaped profile. For 
simplicity we shall treat the case of interference patterns due to just a pair of slits of width d and distance D, even 
though an analytical solution can be obtained also for a grating composed of several slits. 

Let us consider an initial wave function split by two Gaussian slits of standard deviation <j x and centered at 
x = ±D/2 (we can typically choose 6a x = d) 



ipo{x; k x ) 



D 



+ ip[ x 



A k x x 



c 



e + e 



where C is the normalization constant. Inserting Eq. I|32ll in Eq. (|40|l with the initial state (|A1|) we get 

T , \ MC 2 ff + °° , . f iM r 2 , 2 „ , , , n (xq-x'o) 2 
I(x) = J J dx dx' Q exp I— [4- x' 2 + 2x(x' - so)] - \ p °' 



x ex p 



(£0 



i D\ 2 



(Al) 



whence 

7 fx) = . exp 

hT^W) 



D 2 /AfV2Ma x 



Q(T) V hT 



cosh 



/ y/2M<7 x A 2 Dx 
V ) QiT) 



-exp 



L> 2 



2Q{T)P 



MDx 



hTQ(T) 



(A2) 

where Q(T) = 1 + [2ALr 2 /(frT)] 2 + {2a x /£) 2 . 

For T ^> M(T x D/h and £ ^> er^ the term Q(T) — > 1, the first exponential in Eq. (|A2() reduces to 
exp{— [v / 2Mxcr 2 ,/(fi,r)] 2 } , whence x maa; w hT/(Ma x ), and thus the argument of the hyperbolic cosine is close 
to zero. It follows that Eq. I|A2J) is well approximated by (see Eq. JHJ) 
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which coincides with Eq. I|47|l evaluated for N = 2. 



APPENDIX B: APPROXIMATION QUALITY 
TEST FOR FULLERENE EXPERIMENT 

In this appendix we evaluate the error committed 
introducing the zero-order approximation of (|45|) in 
Eq. (|44H . referring to experimental conditions reported 
in 0. 

Since \ (s' — s)D\ < 2£, we will test the previous approx- 
imation in the most unfavorable situation, according to 
exp[(£' — £)/£] ~ 1. To this aim, it is useful to introduce 
the integrals 

r+oo 

1(1, k)= / d^fO^e 1 ^, 



A(k) 



where k = Mx/(hT). Zero-order approximation of 145|) 
can be checked by evaluating the relative displacement 



of the square modulus of the previous integrals 
\I(£M 2 -\A{k)\ 2 



R(£,k) 



mkw 



The less is the value assumed by R(£, k), the better is the 
quality of the zero-order approximation of (|45H . 

By making explicit tp(£) with the characteristic func- 
tion defined in the interval [— d/2, d/2], a straightforward 
calculation leads to 



\l(£,k)\ 2 
\A(k)\ 2 
whence 
R{£,k) 



1 , £2 + p [cosh(d/£) - cos(fcd)], 



d 2 sinc 2 (fc<i/2) 



p- [1 - cos(fcd)], 



1 



{£kf 



1 — cos(fcd) 



cosh(d/£) — cos(kd) 



(Bl) 
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In the limit of small k (i.e., for positions x « 0, close to 
the maximum of intensity) , we get 



(d/£) 2 

lim R(£, k) = 1 - -p — \-f-rh: -r 

k^o V ■ ' 2 cosh (d/i) - 1 



12 



0.011 



This means that in correspondence of the greatest inten- 
sities the approximate expression just moves away from 



the correct one of about 1%. 

Since R{£, k) is an increasing function of k, and thus 
of \x\, then the most unfavorable case takes place at the 
edge of the interference pattern, where, however, the in- 
tensity is negligible. By using fc max ~ d^ 1 in evaluating 
I|B1(I . it results that R(£, k max ) does not significantly dif- 
fer from evaluation (IB2I) . 
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scribed by the reduced density matrix p(r, r') by not- 
ing that: 1) the analysis in Q extends to "conditional 
wave functions"; 2) the reduced density matrix p(r, r') 
is indeed the projector onto the pure state defined by 
the conditional wave function averaged with respect to 
the "quantum equilibrium distribution" . The notions of 
"conditional wave function" and "quantum equilibrium 
distribution" have been introduced in [16| and are cru- 
cial for a proper understanding of the empirical import 
of Bohmian mechanics. The only delicate point is current 
positivity. In order to ensure that the statistics of escape 
time and position is J, current positivity should hold for 
all the conditional wave functions whose average leads to 
P(r,r'). 

Eq. 1201 differs from the analog in y| and |1m| for numer- 
ical constants, which here have been corrected. 
The values of pressure and cross section were kindly sup- 
plied by Dr. Olaf Nairz. 

A recent proposal could lead to localization factors 



smaller by a factor of 2ir 30] . In this case decoherence 
effects, which are already quite less relevant with respect 
the loss of coherence due to beam preparation, should be 
further negligible. 

[56] Eq. holds under the conditions (gfl, T > r with 
r = MAx£/h, and d/i < 1. These conditions imply 
a large superposition of the wave packets on the screen. 
In fact, the ratio 1Z between the separation of the most 
distance slits (~ Ax) and the size of the pattern I(x) 
(~ \L/2nd = hT/Md) becomes TZ < d/l 

[57] Note that condition 1461 1 . £ S> d, is obviously consistent 
with Ij48^ . i > D, being d < D. In particular, for fullerene 
experiment, condition I46L via definition 139H . leads to 

T « 2iff[l"W) 2 ] ~0.3s. 
[58] In other words, we have taken the <^ s (a;)'s in J3TJ to 
be characteristic functions with supports on the effective 
width of the slits. Numerical simulations show stability 
with respect to other equivalent choices for these func- 
tions. 



